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Abstract 

The goal of the present paper is to propose an enhanced ordinary differential equations solver 
by exploitation of the powerful equivalence method of Elie Cartan. This solver returns a target 
equation equivalent to the equation to be solved and the transformation realizing the equiva- 
lence. The target ODE is a member of a dictionary of ODE, that are regarded as well-known, 
or at least well-studied. The dictionary considered in this article are ODE in a book of Kamke. 
The major advantage of our solver is that the equivalence transformation is obtained without 
integrating differential equations. We provide also a theoretical contribution revealing the re- 
lationship between the change of coordinates that maps two differential equations and their 
symmetry pseudo-groups. 
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1. Introduction 



Current symbolic ODE solvers make use of a combination of Lie symmetry methods 
and pattern matching techniques. While pattern matching techniques are used when the 
ODE matches a recognizable pattern (that is, for which a solving meth od is already imple- 
mented) , sy mmetry methods are reserved for the non-classifiable cases (|Cheb-Terrab et al 



1997, 199- 
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Fig. 1. General flowchart of typical ODE-solver. 

Nevertheless, in practice, these solvers often fail to return closed form solutions. This 
is the case for instance of the equation 

,2 1 

(1) 



y 2 



which admits only one 1-parameter symmetry group. Using this information, present 
solvers return a complicated first order ODE and a quadrature which is quite useless for 
practical applications. 

More dramatically, when applied, to the following equation, these solvers output no 
result 

y" - gaV-fig 2 *-! = o. (2) 

x b 

This failure is due to the fact that the above equation does not match any recognizable 
pattern and has zero-dimensional point symmetry (pseudo-)group. Thus neither symme- 
try methods nor classification methods works. The same goes for the equation 



y"y + y' - 4 y e + 12 xy 4 - (4 x + 12 x 2 ) y 2 + 4 x 3 + 4 x 2 - 2 a = 



(3) 



Our solver is designed to handle such differential equations. It returns a target equation 
equivalent to the equation to be solved and the equivalence transformation. The target 
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ODE is a member of a dictionary of ODE, that are regarded as well-known, or at least 
well-studied. The dictionary considered in the article are ODE in a book of Kamke. 
However other ODE could be added to the dictionary without difficulty and this is an 
advantage of the method. 

For the equation (1), we obtain the Rayleigh equation (number 6.72, page 559 in 
Kamkc's book) 

y" + y' 4 + y = 

and the change of coordinates 

(x, y) -> (x, —). 

For the equation (2), we obtain the first Painleve equation (number 6.3, page 542 in 
Kamke's book) 

y" = 6y 2 + x 

and the change of coordinates 

(x, y) -> (-, y)- 

x 

The equation (3) is mapped to the second Painleve equation (number 6.6 in Kamke's 
enumeration) 

y" = 2y 3 + yx + a 

under the transformation 

(x, y) -> (2x, y 2 - x). 
To incorporate such changes of variables, one needs to understand the equivalence 
problem : Given two ODE 

E f : y" = f(x,y,y'), Ef : y" = f(x,y,y') 

and an allowed Lie pseudo-group of transformations T$ acting on the variables x and y (for 
real-life reasons, we restrict ourselves to second order ODE. However, what follows re- 
mains true for any order). We shall say that Ef and Ef are equivalent under the action 
of the pseudo-group r<E> if there exists a change of coordinates ip G T$ that maps one 
equation to another. This will be denoted by 

Ef — (p*(Ef) and tp E T$, (4) 

or in abridged form Ef ~$ Ef. As we shall see, the system (4) is PDE's system. One can 
consider the equivalence of an equation with its self : A symmetry a of the equation Ef is 
a solution of the self-equivalence obtained by setting Ef — Ef in the a PDE's system (4). 
The solutions of this self-equivalence problem form a Lie pseudo-group, the symmetry 
pseudo-group, which will be denoted by SE f ,<s>- 

In practice, one distinguishes two possible situations in the computation of the change 
of coordinates. First, the input equation (the equation to be solved) and the target 
equation are known. This is an online computation. In the second situation, considered 
here in the construction of our solver, only the target equation (an equation from Kamke's 
list) is known and we look for the change of coordinate that maps the generic equation 
to the target one. So, assume that Ef is a generic differential equation and Ef is fixed 
equation that falls within the effective differential algebra i.e. Ef is given by equalities 
between differential polynomials with rational coefficients. 

Crucial in the construction of our solver, is the establishment of the relationship be- 
tween the change of coordinates and symmetry pseudo-groups. In particular, one might 
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ask under which conditions the change of coordinates can be obtained without integrating 
differential equation ? To the authors's knowledge, this is the first time in equivalence 
problem theory that such questions are investigated. The answer, which constitutes the 
theoretical contribution of the paper, can be be summarized as follows (we emphasis on 
the fact that what follows remains true for any order) : 
(i) The number of constants appearing in the change of coordinates (p S r$, map- 
ping Ef to Ej, is exactly the dimension of Se^,®- This implies that when this 
dimension vanishes the change of coordinate can be obtained without integrating 
differential equations. Also, we have dim(S l £; /i $) = dim(S , Sj,$). 
(ii) In the particular case when dim(5 l £; / -,*) = 0, the transformation (p is algebraic 
in / and its partial derivatives. The degree of this transformation if is exactly 
card(5 , £ / - i $). In this case, the symmetry pseudo-groups SE f ,<s> and Set,® have finite 
cardinals. However, they need not to have the same cardinal. 
The simple fact dim(S , £; /i $) = dim.(SEf,$) allows us to construct a powerful hashing 
function which significantly restricts the space of research in kamke's list. For this reason, 
we use 7 possible types of transformations $1, • • • , $7 (see table 1 page 14). We pre- 
calculate to each target equation a signature index, that is, the dimensions of the 7 
symmetry pseudo-groups associated to the 7 types of transformations (this calculation 
is done without any integration). If two differential equations are equivalent then their 
signature indices match. 

The transformation tp in (ii) can be obtained using differential elimination algorithms. 
This is explained in third section. Unfortunately, such approach is rarely effective due 
to expressions swell. In order to avoid this, we propose in section 4 a new method to 
pre-compute the transformation ip in terms of differential invariants (we do this for each 
target equation Ej in Kamke's list). These invariants are provided by Cartan's method. 



2. Equivalence problem 

The equivalence problem is the study of the action of a given pseudo-group of trans- 
formations on a given class of differential systems. In the algebraic framework, this action 
is viewed as the action of a X>-groupoid $ on a diffiety £. 

2.1. Equivalence problems and groupoids 

Recall that a diffiety (see A.l) is the set of formal Taylor series which are regular 
solutions of a finite PDE's system. It is a pro-algebraic variety, fibered over an algebraic 
variety X and which will be denoted by ir : £ — > X. The projection of a Taylor se- 
ries £ £, of a function / : X — > U, is the expansion point x G X . The coordinate 
ring of a diffiety is a reduced finitely generated differential algebr a. The autom orphisms 
group of the diffiety is the set of the contact transformations (see Olver ( 1993f) ) from £ 
to £. 

A P-groupoid $ is a diffiety formed by invertible Taylor series and closed under the 
composition (see section A. 4). A X>-groupoid acting on a manifold X, is a subset of the 
space of infinite invertible jets J£°(A, X). The Taylor series of contact transformations 
of a diffiety £ form a X>-groupoid that acts on £ and which will be denoted by aut(£). 



4 



Given a diffiety £ fibred over X and a X>-groupoid <£> acting on X. An equivalence 
problem is the action of $ on the diffiety £ , that is, an injective representation 1 1, i.e. an 
injective morphism of P-groupoids 

p : $ — > aut(£). 

Example 1 (2nd order ODE, x = x + C,y = rj(x,y)). The infinite dimensional T>- 
groupoid $ := $3 acts on the points (x,y) € J°(C, C). The corresponding Lie defining 
equations are (we set x = t;(x,y)) 

6 = 1, Hy = 0, Vy + 0. (5) 

The action of $ on the jets space J°(C,C) is prolonged (see appendix A. 5) to an action 
on the first order jets space X := J 1 (C,C). In the coordinates (x,y,p = y'), this action 
reads 

x = £, y = *n, p = f] x +pr] y . (6) 
The equivalence condition (4) is obtained by prolonging the action of $ on the second 
order jets space J 2 (C, C) and by setting y" = f(x, y,p). Thus, we obtain 

£r = 1, £,y = 0, £,p = 0, W T] y = 1, 

x = £, y = *n, p = Vx+PVy, (7) 

f{x,V,p) = Vxx + 2pr)xy + P 2 Vyy + f(x,V,P) Vy 

This action of $ on the equation Ef is viewed as an action on the Taylor series J^f, 
in other words, as an action on the trivial diffiety £ :— 3 oc (X, C) (fibred over the mani- 
fold X). The coordinates ring of X is C[X] := C[x, y,p] and the coordinate ring of £ is 
the ring of differential polynomials 

C[£] :=C[,, y) p] { /}withA = {^|,|}. 



3. Differential-algebraic approach 

The aim of this section is to use differential elimination to solve the equivalence prob- 
lem (4) when the target function f(x,y,p) is a Q-rational function, explicitly known 
and the symmetry pseudo-group of Ej is zero-dimensional. The reader can find in the 
appendix A.l a brief introduction to differential algebra. 

3.1. The self- equivalence problem 

The system (4) is fundamental and can be treated by two different approaches : brute- 
force method based on differential algebra (this section) and geometric approach relying 
on Cartan's theory of exterior differential systems (the next section). While it is classi- 
cally known that the existence of at least one transformation (p : X — > X can be checked 
by computing the integrability conditions of the system (4), whic h is completely algo - 
rithmic whenever the functions /, / : X — > C are explicitly known ( Boulier et all [l995h . 



there is no general algorithm for computing closed form of tp. We shall show that if 
the symmetry pseudo-group of Ei is zero-dimensional, the transformation ip is obtained 
without integrating any differential equations. 



Elie Cartan call prolongement holohedrique such injective representation. 
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Definition 1 (Symmetry pseudo-group). To any differential equation Ei and any T>- 
groupoid <3? that acts on (x, y), we associate the £>-groupoid Se^.® formed by the formal 
Taylor series solutions of the self-equivalence problem 

E f = o*{E f ) and a e F$. (8) 

The symmetry pseudo-group SW# := YSej,^ is the set of C°°-functions a : X — > X 
that are local solutions of the Lie defining equations (8). 

Example 2. Consider the £>-groupoid < fr := $3 f rom e xample 1 and the Emden- Fowler 
equation Ej (number 6. 11, page 544 in (<KamkeL[l94l ) 

y" = — 2 - (9) 

xy z 

The Lie defining equations of the symmetry pseudo-group of Ej are obtained by setting 

f(x,y,p) = — tt and f(x,y,p) — ^ttt in the equations (7). The characteristic set (see 
xy 1 xy 1 

definition 8 page 16) of these equations is 

c a = U = v-, t = y\ x = x\. (10) 



y 

This PDE's system is particular : it contains only non-differential equations. We have 
dim C CT = and degC CT = 3. The symmetry pseudo-group Sjr := Sej,<s> is actually a 
group with 3 elements: 

Sf = {{x,y,p) -► (x,Xy,Xp) A 3 = 1} 

3.2. Equivalence problem with fixed target 

Assume that f(x,y,p) is a fixed (Q-rational function. 
Example 3. Consider again the equivalence problem of example 1 : 

p-y~x- py y = 0, 

y xx + 2py xy + p 2 y yy + f y y - f(x, y, p) = 0, (11) 
X x — 1 = 0, x y = 0, x p = 0, y p = 0, y y ^ 0. 

These equations constitute a quasi-linear characteristic set w.r.t. the elimination ranking 

e/ >- ep >- eg y Qx. 

Consequently, the associated differential ideal is prime. This fact can be generalized to 
any P-groupoid $ defined by quasi-linear characteristic set (see proposition 8). 

Proposition 1. The PDE's system (4) (where f(x,y,p) is a fixed rational function) is 
a quasi-linear characteristic set C w.r.t the ranking 6/ >- Op >- 6{y, x}. 

3.3. Brute-force method 

Using Rosenfeld-Grobner we compute a new characteristic set C of the PDE's 
system (4) w.r.t. the new elimination ranking Q{p, y, x} >- 0{/}, which eliminates the 
indeterminates {p, y, x}. We make the partition of C as in the formula (A. 4) 

C = C f UC v (12) 
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where C s := C n ®[x,y,p]{f} and C v :=C\C f . 

Proposition 2. The transformation ip does exist for almost any function / satisfying the 
PDE's system associated to the characteristic set C/. The function x = ip(x) is solution 
of the PDE's system associated to C v . 

Definition 2. When dimC^ = 0, the algebraic system associated to C v is called the 

necessary form of the change of coordinates x = y(x). 

Example 4. Consider the equivalence problem of example 1. Suppose that the target Ej 
is the Airy equation 

y" = xy. 

In this case, Rosenfeld-Grobner returns C v and C/ resp. given by (13) and (14) 



Vxx = -fyy+pfpVy - \v 2 f V pVy + yf v -\ vfx P - \yfppf + \yfy - \ypfyp 

Vxy ^2'fpVy 2Pfp>pVy 

Vyy = -\fppVyi ( 13 ) 
Vp = 0, 

x = fy ~ 2~fxp ~ 2~fppf 4-fp ~ 2.PfvP 



fxxp 2 f X y + fpfxp 2 + p fyyp fppfx ~\~ ' ' ' 

fxyp = 2 fyy — pfyyp — fyppf — fppfy + fpfyp 
fxpp = fyp ~ Pfypp 

fppp = o. 

We have dim C v = 3 which means that the transformation x = <^(x), when / satisfies C/, 
depends on 3 arbitrary constants. 

Example 5. Assume now that the target equation Ej is 

y" = f- 

Rosenfeld-Grobner returns C v and Cj resp. given by (15) and (16) 



JJ 2 — 12 (^./y % f xp 2 fppf 2pfyp + fp ), 

(15) 
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fxxxp — (4 fy 2 f x p 2 fppf 2 pfyp ~\~ fp ) X 

(24p 2 f yp fy — 2Ap 2 fyy.fp.fyp + ' ' ' + 12p 2 fx.yp + 12 / fyp — 8 ,/pp / ) 
fxxyp = (4 /j/ — 2 / xp — 2 fppf — 2 p/yp + fp ) x 

{~^P.fyppf.fpfyp ~ 4 /;jp /yp + • • • + Qpfpfy,ypfppf + 2p 3 /yj / p ) 

2 ' ' 

fxyyp (4 /y — 2 / X p — 2 fppf — 2 pfyp + /p ) X 

(2 ,/j/pp / — 2 fppfyppfpfyp + ' ' ' + 4 /j,p / X p + 4 /yp p + 16 fyppfy ) 

/xpp — /j/p — P.fypp 

fppp = 0. 

Consequently dimC^ = and degC^ = 2. Thus, 95 is an algebraic transformation of 
degree 2, given by equations (15). 

5.^. From equivalence problem with determined target Ej to the self- equivalence problem 

Consider the characteristic set C — Cf U Cp associated to the equivalence problem 
with determined target Ej defined by (12) and computed w.r.t the elimination ranking 
&{p, y, x} y 8{.f }. On the other hand, consider the characteristic set Co- associated to 
the self-equivalence problem (8), computed w.r.t the orderly ranking on {p, y, x}. By 
definition, we have 

C C Q(x,y,p){x, y, p, /} 
C a C Q(x,y,p){x, y, p} 

We obtain C a from C by requiring that the two functions / and / arc equal. The set C 
is specialized by substituting the symbol / by the value f(x, y, p). After specialization, 
the differential system Cf constraining the function / is automatically satisfied since 
there exists at least one solution x = er(x) of the problem, namely a = Id. 

Lemma 1. The two characteristic sets C v and Co- have the same dimension and the 
same degree in the zero-dimensional case. 

Proof. For each equation Ef equivalent to the target equation Ef, denote by <!>y j the 
difhety defined by C v . The £>-groupoid Sej,<s> acts simply transitively (see Figure 2) on 
the diffiety &fj, i.e. 

Vyo, <P € T$fj, 3la e TSej,<s>, ip = a o tp . 

Jx/ — — - Jx-af 




39. f 

Fig. 2. Simply transitive action of Sej,<s> on <&f j where x = <£o(x) and x = ip(x) 



8 



Every tpo G TQjj, define a bijective correspondence 5b /i $ — > 

J> — £V = ° C??Vo), (a G S Ef ,*). 

In fact, according to the Taylor series composition formulae, the one-to-one correspon- 
dence between the two algebraic varieties 5b, t $ and $y j is bi-rational. Consequently, 
these two varieties have the same dimension and the same degree in the zero-dimensional 
case. The same goes for the two characteristic sets C v and C a defining these varieties. □ 

We are now in position to announce the main theorem of the paper (recall that / is 
assumed Q-rational function of its arguments and by definition, the degree of an algebraic 
transformation x = y>(x) is the generic number of points x when x = (x,y,p) G C 3 is 
determined) : 

Theorem 3. The following conditions are equivalent 

(1) &m{C 9 ) = 0, 

(2) dim(C CT ) = 0, 

(3) <tim(S Ef ,<s>) = 0. 

(4) caid(SEf,<s>) < oo. 

In this case, card 5b = deg(C v ) = deg<^. 

Remark 1. When the transformation x := y(x) is locally bijective, but not glob- 
ally, S , £; /j $ and 5bj,$ need not to have the same degree. Indeed, consider again the 
groupoid $3 and the equations 

„ 6y 4 + x- 2y' 2 _„ 2 _ 
y = and y = 6y + x 

which arc equivalent under (x — x, y = y 2 ). The corresponding symmetry group are 
respectively given by 

Sej,® = {(x,y) -> (x, Ay) | A 2 = 1} and S Ef -,<i> = {Id}. 
They have the same dimension but different cardinality. 

3.5. Expression swell 

In practice, the above brute-force method, which consists of applying Rosenfeld- 
Grobner to the PDE's system f 4), is rarely effective due to expressions swell. Most of 



the examples treated here and in (|Dridil . 12007). using our algorithm ChgtCoords, cannot 
be treated with this approach. 

It seems that the problem lies in the fact that we can not separate the computation 
of C v from that of Cf which contains, very often, long expressions (observe that, since C v 
is computed in terms of / then there is no need to compute the equivalence conditions Cf). 

An other disadvantage of the above method is that we have to restart computation 
from the very beginning if the target equation is changed. In the next section, we propose 
our algorithm ChgtCoords to compute the transformation <p alone and in terms of differ- 
ential invariants. These invariants are provided by Cartan method for a generic / which 
means that we have not to re-apply Cartan method if the target equation is changed 
and a big part of calculations is generic. Furthermore, the computation of <p in terms of 
differential invariants significantly reduces the size of the expressions. 
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4. Cartan's method based approach 

We refer the reader to (jCartanl[l953l ; lHsu and Kamranl . [l^^lOlveii[l995l ; lNeutl . [2005 



Dridil , 120071 ) for an expanded tutorial presentation on Cartan's equivalence method and 



application to second order ODE. 

When applied, Cartan's method furnishes a finite set of fundamental invariants and a 
certain number of invariant derivations generating the field of the differential invariants. 

Example 6. Consider the equivalence problem of example 1. The PDE's system (4) 
reads (p — y') 



dp - f(x,y,p)dx^ 


U 


«2 


o\ 


' dp - f(x,y,p)dx^ 


dy — pdx 


= 







dy — pdx 


{ dx ) 


v° 





V 


{ dx ) 



uij S(a) Uf 

with det(S'(a)) ^ 0. In accordance with Cartan, this system is lifted to the new linear 
Pfaffian system 

S(&) loj = S(a) LUf. 

After two normalizations and one prolongation, Cartan's method yields (an e-structure 
with) three fundamental invariants 



1 



h = - 7 (f P r-f y +M 1 



4 wf/ Jy 

and the invariant derivations 

ld_ 
a dp' 
ld_ 
a dy 



h = 



x 2 = 



X 3 = D a 



2 a dp 2 fpp 8a 



2a 2 
d 





2 fpa da> 



h = 



X A 



fyp &x fp 

2a 



da 



(17) 



(18) 



where a = a% and as usual D x — ^ + f( x i ViP)~§j; denotes the Cartan vector field. 

When dim(5£; / - i $) = 0, the additional parameters a and a can be (post)normalized by 
fixing some invariant to some suitable value. 

Proposition 3 (jOlven (|l995l )U The symmetry groupoid is zero-dimensional if 

and only if there exist exactly three functionally independent specialized invariants 
h[f},h[f},h[f}:X^C. 

Note that the invariants I\ [/],•■ ■ , ^3 [/] are functionally independent if and only if 
dl\ [/] A ■ ■ • A d/3 [/] 7^ and if the function / is rational, then the specialized invariants 
/[/] : X — > C are algebraic functions in {x,y,p). 

4-1. Computation of ip 

Suppose that the symmetry groupoid 5^-,$ is zero-dimensional. Then, according to 
theorem 3, there exists 3 invariants Fj. := /&[/], 1 < k < 3 such that the algebraic (non 
differential) system 

{Fi(x) = Ji,F 2 (x) = h,F 3 {5t) = I 3 } (19) 
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is locally invertible and has a finite number of solutions 

x = F- 1 (/ 1 ,...,/ 3 ). 
The specialization of I\ , . . . , I3 on the source function / yields 

x = F-\l 1 [f],...,I 3 [f}). 



(20) 



(21) 

The main idea here is that the inversion (20) is done by computing a (non differential) 
characteristic set C for the system (19) w.r.t. the ranking {x, y, p} y {Ii, I2, J3}. Now, 
the most simple situation happens when deg(C) = 1. In this case, the necessary form of 
the change of coordinates <p is the rational transformation defined by C . 

Example 7. Consider the equi valence problem of the example 1 and the target equa- 
tion Ej introduced by G. Reid (jReid et all Il993l ) 

4f 



y" = y -r 

x 



The following invariants are functionally independent (where Ii-j-.-k — Xk ■ ■ ■ Xj(ij) ) 



h = — + 8 — 

Ax 2 x 3 ' 



'1:3 



-3 x — 48 y + 16 px 
2X 1 ' 



'1;23 



-20- 



ax^ 



'l;31 



ax° 



We normalize the parameter a by setting I123 = —20. The characteristic set C is 



3 

32 



512 



h 



31 



4096 



:^l;3^1;3l'" 



{ V = ^1-31 

^ J 256 ' 



X = gll;31, 



——I 3 

4096 /l;3i: 



which gives the sought necessary form of ip. As a byproduct we deduce that the symmetry 
group Se t ,<s> = {Id}. 

When deg(C) is strictly bigger than 1, we have two cases. First, deg(C) = deg(SE f ,<s>) 
and then <p is the algebraic transformation defined by C. Second, deg(C) > deg(5 l s / -,<i>)- 
In this case, to obtain the transformation tp, we have to look for 3 other functionally inde- 
pendent invariants such that the new characteristic set C has degree equal to deg(S l £; / -,<i>). 

Example 8. Cons ider the equi valence problem of example 1 and the target equa- 
tion number 6.9 in (|Kamkelll94i : 



y 

The corresponding symmetry group is 



y +xy. 



Sej,9 = {{x,y) -> (x, Xy)\X 2 = 1} . 

One can verify that I\, Ii ; i3 and Ii ; i33, when specialized on the considered equation, are 
functionally independent. The associated characteristic set C is 
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(Ax 2 + 2hx - 3/i;33 - 2/i 2 ) 



3(/i ;3 + 1) 



3 T o 3 , 3 _ 3 T n 3 , 1 o 3 

-hx 2 + -I 1;33 x - -J 1;3 - -7 1;3 2 + -/i i33 /i + - g> 



which has degree 6, strictly bigger than the degree of the symmetry group. However, if 
we consider the invariants K\ := /i ; 23 3 /7i ;3 i, K 2 := h-234,/ h-,si and K 3 := /i : 23i/7 2 ;3 i, 
we obtain 



p = -Kiy, 
1 

S 

x = --K z + Kx. 
6 



y 2 = 7^3 



That is, the necessary form of <p since this new set has degree two. 
4-2. Heuristic of degree reduction 

We have to remark here that searching invariants giving the required degree (as in 
example above) is not an easy task, although algorithmic. This is simply because the 
algebra of invariants can be very large. For this reason we provide an important heuristic 
which enables us to obtain the desired degree even for a "bad" choice of invariants. 

Example 9. Consider the Emden-Fowler equation (9) and the X>-groupoid of trans- 
formations $ 3 . We have already computed the corresponding symmetry groupoid. The 
specialization of the invariants I\, Ji ; i3 and 7i ; i33 gives three functionally independent 
functions. As explained above, we obtain the following characteristic set computed w.r.t. 
the ranking p>-y>~x>~Ii>~ Ji ;3 >- 7i ;33 

/3 17 1;33 , l/i ;3 2 \ __ l/i i3 _ 

"- (J --^ + l¥)- x -^ (22) 



9 Ji 6 - 8/ 1;3 z + 67 1;33 /i J MS - 8I h ^ + 67 1;33 /i 

Comparing with the 2?-groupoid of symmetries (10) we deduce that, in contrast to y, 
the degree of x must be reduced to one. This can be done in the following manner. First, 
observe that the Lie defining equations of $3, more exactly x p = 0, implies that X\ (x) = 
where X\ — is the invariant derivation (18). Now, differentiate the last equation of 
the characteristic set, which we write as x 2 = Ax + B, w.r.t the derivation X\. We find 
A.^x + B : i = 0. The coefficient of x in this equation, which is invariant, could not vanish 
(since it is not identically zero when specializing on the Emden-Fowler equation). We 
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£>l 

obtain x — — -p- or explicitly 



;1 



2 Kh^+hKg 



with K = 



(23) 



9/j 3 -8I 1;3 2 + 6I 1;33 I 1 



The necessary form of the change of coordinates <p is then given by (23) and the first two 
equations of (22). 

The above reasoning can be summarized as follows 
Procedure ChgtCoords 



Output : x = iys(x) the necessary form of the change of coordinates 

1- Find 3 functionally independent invariants -M/], -^[/D defined on X. 

2- Compute a char, set C of the algebraic system (20). 

3- If dcg(C) = 1 then Return C. 

4- Compute Sej,& with Rosenfeld-Grobner. 

5- WHILE deg(C) ^ deg(5 B/ -,$) DO 

Reduce the degree of C. 
END DO 

6- Return C. 



5. The solver 

5.1. Pre- calculation of ip 

5.1.1. The first step : the adapted V-groupoid 

Let $i, . . . , $7 C JJ°(C 2 , C 2 ) denote the D-groupoids defined in the table 1 page 14. 
It is not difficult to see that $i C $ 3 C $ 5 and $ 2 C $4 C $ 6 and finally $ 5 , $ 6 C $7. 

Definition 4 (Signature index). The signature index of Ef is 

sign(£/) := ((di, d 3 , d 5 ), (d 2 , d 4 , d 6 ), d 7 ) where d 4 := dim5 E/ ,* i , 1 < i < 7. 

Clearly, (di < d3 < ds < d7) and (d 2 < d 4 < d6 < d7). Recall that the calculation of 
theses dimensions does not require solving differential equations. 

Definition 5. We shall say that the signature index sign(_E/) matches the signature 
index sign(Sj) if and only if 



Input : Ef and <f> such that dim(5 , £^ i $) = 



d7 = d7 and (si = S\ or s 2 = s 2 ) 
where s\ and s 2 stand for (di,d3,ds) and (d 2 ,d4,dg) resp. 
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Transformations 


Equation number according to Kamke's book 


$1 


x = x, y = ??(a;, y) 


1, 2, 4, 7, 10, 21, 23, 24, 30, 31, 32, 40, 42, 43, 45, 47, 50 




x = x + C, y — rj(x, y) 


11, 78, 79, 87, 90, 91, 92, 94, 97, 98, 105, 106, 156, 172 


$5 


x = £(x), y = 17(1,1/) 


Null 


$ 2 


S = €(x,y), y = y 


81, 89, 133, 134, 135, 237 


$4 


x = £(x,y), y = y + C 


11, 79, 87, 90, 92, 93, 94, 97, 98, 99, 105, 106, 172, 178 


$6 


x = i{x,y), y = v(y) 


80, 86, 156, 219, 233 


$7 


x = £(x,y), y = rj(x,y) 


3, 5, 6, 8, 9, 27, 44, 52, 85, 95, 108, 142, 144, 145, 147, 
171, 211, 212, 238 



Table 1. Adapted groupoids for certain equations from Kamke's list 



Definition 6. Two second order ODE Ef and Ej are said to be strongly equivalent if 
3$ G {$1, • ■ • ,$7}, lip G T$, ip*E f = Ej, dim <SW $ = 0. 

Lemma 2. If Ef and Ej are strongly equivalent then their signature indices match. 

Definition 7 (Adapted 2?-groupoid). A £>-groupoid $ is said to be adapted to the 
differential equation Ej if dim(S l £; / .$) = and $ is maximal among $1, ■ • • , $7 satisfying 
this property. 

The table 1 associates to each equation in the third columr[H] its adapted groupoids. 
For instance, the first Painleve equation (number 3) appears in the last row which means 
that its adapted P-groupoid is the point transformations £>-groupoid $7 . To the Emden- 
Fowler equation, number 11, we associate the Z?-groupoids $3 and $4. In the case of 
homogeneous linear second order ODE (e.g. Airy equation, Bessel equation, Gaufi hyper- 
geometric equation) we prove that, generically, the adapted P-groupoid is $4. 

5.1.2. The second step 

Once the list of adapted D-groupoids $ is known, we proceed by computing the neces- 
sary form of the change of coordinates <p € T$ using ChgtCoords. Doing so, we construct 
a Maple table indexed by Kamke's book equations and where entries corresponding to 
the index Ej are: 

1- the signature index of Ej, 

2- the list of the adapted P-groupoids $ of Ej, 

3- the necessary form of the change of coordinates ip G r<f>. 

For instance, the entries associated to Rayleigh equation y" + y' 4 + y = are: 

1- the signature index ((0, 1, 1), (1, 1, 1), 1), 

2- the 2?-groupoid $1, 

3- the necessary form of the change of coordinates 



2 A more complete list of equations is available upon request. 
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p 



-36 



h-i 



72 + 72 Ji + J 2i i 




< 



?7 :i 



V 



with the normalization I2/I2.1 — 1- Invariants here are those generated by (17) and (18) 
plus the essential invariant x = x. 

5.2. Algorithmic scheme of the solver 

To integrate a differential equation Ef our solver proceeds as follows 

Procedure Newdsolve 
Input : Ef 

Output : An equation Ef in Kamke's list and the transformation ip such that 



(i) Specialize, on Ef, the necessary form of the change of coordinates associated 
to Ef. We obtain tp. 

(ii) If ip € T$ and ip*(Ef) — Ef then return (Ef, ip). 
END DO. 



5.3. Features of the solver 

It is worth noticing that the time required to perform steps (i)- (ii) is very small. In 
fact, it is about one hundredth of a second using Pentium(4) with 256 Mo. Experimented 
on many examples, the total time needed to solve a given equation does not exceed few 
seconds in the worse situations. 

The second feature of our solver is, contrarily to the symmetry methods, neither 
the table construction nor the algorithm of the solver involves integration of differential 
equations. Indeed, even the computation of signature indices is performed without solving 
the Lie equations. 



<p.(E f ) = E f 



1- Compute the signature index of Ef. 

2- Select from the table the list of equations Ef such that sign(Ef) matches sign(£/). 

3- FOR each equation Ef in the selected list DO 
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A. Appendixes 



A.l. Differential algebra 



The reader is assumed to be familiar with the basic notions and notations of diffcr- 
ential algebra. Reference books are dmtl Il950h a nd (|KolchiiiL [l973h . We also refer to 
(|Boulier et all Il995t iHubertl . l2000l ; iBoulieii l2006h . Let U = {m, ■ ■ ■ , u n } be a set of 
differential indeterminates. k is a differential field of characteristic zero endowed with the 
set of derivations A = {di, • ■ • , d p }. The monoid of derivations 



9 := {d^d. 



C*2 

2 



■d a " 



Oi\, 







(A.l) 



acts freely on the alphabet U and defines a new (infinite) alphabet QU. The differential 
ring of the polynomials built over QU with coefficients in k is denoted R = k{£/}. 
Fix an admissible ranking over QU. For / G R, ld(/) G QU denotes the leader (main 
variable), If G R denotes the initial of / and Sf £ R denotes the separant of /. Recall 
that Sf — ^ where v = ld(/). Let C C R be a finite set of differential polynomials. 
Denote by [C] the differential ideal generated by C and by \J\C\ the radical of [C]. Let 
Hp '■= \If I / € C}L){Sf I / G C}. As usual, fulLrem is the Ritt full reduction algorithm 



(|KolchirJ . Tl973h . If r =full.rem (/, C) then 3h G i?^ 3 , /i/ = r mod [C]. Then the reduced 
form is defined by reduced_form(/) := r/h. 

Definition 8 (Characteristic set). The set C C R is said to be a characteristic set of 
the differential ideal c := \/[C] : H^ if 

(1) C is auto-reduced, 

(2) / G c if and only if fulLrem (/, C) = 0. 



Definition 9 (Quasi-linear characteristic set). The characteristic set C C R is said to 
be quasi-linear if for each / G C we have deg(/, v) — 1 where v is the leader of /. 

Proposition 4. When the characteristic set C is quasi-linear, the differential ideal 
c := '■ He C i? is prime. 

A. 2. Taylor series solutions space 

Let k := C(xi, ■ • ■ ,x p ) be the differential field of coefficients endowed with the set 
of derivations | , • • ■ , j • Let C be a characteristic set of a prime differential ideal 
c c R. We associate to C the system 

(C = 0,F c ^0) (A.2) 

of equations / = 0, / G C and inequations h 7^ 0, /i G -He- 

Definition 10 (Taylor series solution). A Taylor series solution (with coefficients in C) 
of the PDE's system (A.2) is a morphism /x : R —> C of (non differential) C-algebras such 
that 

[C] C ker n and ff c n ker \i = 0. 
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The source of the Taylor solution fj, is s(/x) := (^(xi), . . . , ^{x p )) G C p and the target 
is t(ji) := (/i(tii), . . . ,fj,(u n )) G C™. The difhety associated to the characteristic set C is 
the set of the formal Taylor solutions of the system (A. 2). 

The dimension of the solutions space of (A. 2) is the number of arbitrary constants 
appearing in the Taylor series solutions fi when the source point x := s(/j) g C p is 
determined. Let K be the fractions field Frac(i?/c). Recall that the transcendence de- 
gree of a field extension K/k is the greatest number of elements in K which are k- 
algebraically independent. The degree [K : k] is the dimension of K as a k- vector space. 
When tr deg(if/k) = 0, the field K is algebraic over k and [K : k] < oo. If / G C, we 
denote rank(/) :— (v,d) where v := Id / and d := dcg(f,v). Let 

rankC:={rank(/) / G C} 

ldC:={ld(/) l/GC} 
dimC := card (QU \ 6 (Id C)) 

degC:= J]deg(/,ld/). 

fee 

Proposition 5. dimC = tr dcg(if/k) is the dimension of the solutions space of (A. 2). If 
dim C = then the cardinal of the solutions space is finite and equal to deg C = [K : k] . 

A. 3. Differential elimination 

Let U — U\ U U2 be a partition of the alphabet U. A ranking which eliminates the 
indeterminates of U2 is such that 

v«i e eu u vt> 2 e eu 2 , v 2 y Vl . (a.3) 

Assume that C is a characteristic set of the prime differential ideal c = \/[C] : 
w.r.t. the elimination ranking QU2 >~ QU\. Let R\ := k{U{\ be the differential polynomi- 
als fc-algebra generated by the set U\. Consider the set Ci := C fl R\ and the differential 
ideal ci := c fl R\. 

Proposition 6. C\ is a characteristic set of Ci. 

Consider the differential field of fractions K :— Frac(i?/c) and denote by a : R — > K 
the canonical fc-algebra morphism. Let K\ be the differential subfield of K generated 
by the set a{R\). Then K\ is the fraction field associated to the prime differential ideal 
Ci := c fl R\. The partition of the characteristic set 

C = dUC 2 (i.e. C 2 :=C\Ci). (A.4) 

enables us to study the field extension K/K\. 

Proposition 7. tr deg(if/ifi) = dimC 2 . If dimC 2 = then [K : K x ] = degC 2 . 
A.4- Groupoids 

Definition 11 (Groupoid). A groupoid is a category in which every arrow is invertible. 
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Let (3>,X, o, s, t) be a category. Each arrow £ $ admits a source s(ip) £ X and a 
target t{f) 6 X which are objects of this category. The composition 952 °f the two 
arrows p\ and ip? are defined when t(y>i) = s (<^2)- 

If <f> is a groupoid, for each arrow ip G <I>, there exists a unique inverse arrow p~ x such 
that i/j^ 1 o ip = Id s ( v ) and o (,c _:L = Id t ( ¥ ,). 

Let X and [7 be two manifolds and x E X. The Taylor series up to order q (i.e. the 
jet of order q) of a function / : X — > U, of class C 9 , is denoted The Taylor series 
of / about x is denoted j x f or fff. We shall say that x£ X is the source and /(x) 6 17 
is the target of the g-jet 

Example 10. For instance, when X = U = C, we have 

This jet is said to be invertible if f'(x) ^ 0. The jet of the function Id about the point x 
is (x, x, 1, 0, . . . , 0). 

For each integer q 6 N and each x 6 X, we set J|(X, £7) := U/Jx/- We denote 
by J 9 (X, (7) := LLex J x (^> U) the jets space up to order q. We denote by Jl(X, X) the 
submanifold of J q (X, X) form ed by the invertible jets. Recall that Ji(X, X) is a groupoid 
(jOlver and Pohianpelt o, 2006) for the composition of Taylor series up to order q according 



to 



j Ug°f) = (j} {x) g) °(fJ)- (A.5) 

By definition, a 23-groupoid ( Malgranga l200lh $ C 3f(X,X) is a sub-groupoid of 



JJ°(Jf, X) formed by the Taylor series solutions (see def. 10) of an algebraic PDE's system 
called the Lie defining equations. This system contains an inequation which expresses the 
invertibility of the jets. 

The set of C°°-functions ip : X — > X that are local solutions of the Lie defining 
equations of $ is a pseudo-group denoted by T<i>. We define dimr^ = dim<I>/X := dimC 
and, if dimC = 0, cardT$ = card$/X := degC where C is a characteristic set (see 
sect. 3) of the Lie defining equations of $. 

A.5. Prolongation algorithm 

Our aim, here, is to give an efficient way to prolong the action of <J>(°) on J°(C, C) on 
the jets space J™(C, C). For each integer q > 0, define the differential field 

k<«> :=Q{x, y, m,...,y q ) 

and the ring of differential polynomials 

:=k^{x, y, yi,...,y q } 

The differential field k^ is the coefficients field of i£W endowed with the set of derivations 
{cS' J^, . . . , ^-|. Let us assume that the Lie defining equations of are given by a 

quasi-linear characteristic set C R^°\ The P-groupoid acting on J 9 (C,C) and 
prolonging the action of $w is ch aracterized by a characteristic set C . The 
prolongation formulae (jOlver . 19931 ) of the point transformation (x, y) — > (£(x, y), rj(x, y)) 



are of the form 

Vq = Vq{ x > V,---, Vq), 
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where y — r](x, y) if q = 0. The computation of the characteristic set C^ q ' is done 
incrementally using the infinite Cartan field D x := + yi-§- + 2/2^" + • • ■ 

f] q := D x r) q -i ■ {D x £)~ 
c ii) .- u |y 9 _ reduced-form (p 9 , C^ 1 ') } 

Proposition 8. If is a quasi-linear characteristic set of then is a quasi- 
linear characteristic set of w.r.t. the elimination ranking 

Qy q y Qy q -! > y ®{y,x}. 

The previous proposition gives an efficient method to prolong a P-groupoid $ without 
explicit knowledge of (the form of) the transformations tp 6 T$. 
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